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ABSTRACT 

We report the discovery of a transiting planet with an orbital period of 3.05 d orbiting the 
star TYC 7247-587-1. The star, WASP-41, is a moderately bright G8V star (V=11.6) with a 
metallicity close to solar ([Fe/H]= — 0.08 ± 0.09). The star shows evidence of moderate chro- 
mospheric activity, both from emission in the cores of the Call H and K lines and photometric 
variability with a period of 18.4 d and an amplitude of about 1%. We use a new method to show 
quantitatively that this periodic signal has a low false alarm probability. The rotation period 
of the star implies a gyrochronological age for WASP-41 of 1.8 Gyr with an error of about 15%. 
We have used a combined analysis of the available photometric and spectroscopic data to derive 
the mass and radius of the planet (0.92 ± 0.06Mjup, 1.20 ± 0.06i?jup). Further observations of 
WASP-41 can be used to explore the connections between the properties of hot Jupiter planets 
and the level of chromospheric activity in their host stars. 



Subject headings: Extrasolar planets 

1. Introduction 

There is continued interest in finding bright 
stars that host transiting exoplanets because they 
can be accurately characterized and studied in 
some detail, e.g., the mass and radius of the planet 
can be accurately measured. This gives us the op- 
portunity to explore the relationships between the 
properties of the planet and its host star, e.g., the 
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orbital eccentricity, the composition and spectral 
type of the star, the age of the system, etc. Given 
the wide variety of transiting planets being dis- 
covered and the large number of parameters that 
characterise them, statistical studies will require 
a large sample of systems to identify and quan- 
tify the relationships between these parameters. 
These relationships can be used to test models of 
the formation, structure and evolution of short pe- 
riod exoplanets. 

Here we report the discovery by the WASP sur- 
vey of a planetary mass companion to the star 
TYC 7247-587-1. We find that the star is a 
G8V star showing moderate chromospheric ac- 
tivity. The planet, WASP-41 b, is a typical hot 
Jupiter planet with an orbital period of 3.05 d. 

2. Observations 



The WA SP survey i s described in lPoUacco et al 



(|2006l ) and I Wilson et al . (2008) while a discussion 
of our candidate selection methods can be found 



Collier Cameron et al.l ()2007() . iPoUacco et al 
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Fig. 1.— WASP-South photometry of WASP-41. 
Raw data are plotted using small points, phase- 
binned data are plotted using filled circles. Upper 
panel: all data plotted as a function of the orbital 
phase with period P = 3.0524 d. Lower panel: 
data from 2008 plotted as a function of the rota- 
tion phase with period Prot = 18.4 d. 
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Fig. 2. — Faulkes Telescope South z-band pho- 
tometry of WASP-41 (points) with the model fit 
described in Section [5] (sohd line) . 



(j2008l ). and references therein. 

The star TYC 7247-587-1 (WASP-41, ISWASP 
J124228. 50-303823.5) was observed 6767 times by 
one camera on the WASP-South instrument from 
2007 January 20 to 2007 June 22. A further 5637 
observations were obtained with the same camera 
from 2008 January 17 to 2008 May 28. 

The WASP-South lightcurves of WASP-41 
show transit-like features with a depth of approx- 
imately 0.02 magnitudes recurring with a 3.05-d 
period (Fig. [1]). These were independently de- 
tected in the WASP-South photometry from the 
two seasons using the de-tre nding and transit de- 
tection methods described in ICollier Cameron et al 
(j2007l ) , which was taken as good evidence that the 
periodic transit signal was real. The spectral type 
of the star was estimated to be approximately G8 
based on the catalogue photometry available for 
this star at the time. The duration and depth of 
the transit is consistent with the hypothesis that 
it is due to the transit of a planet-like companion 
to a main-sequence G8 star and the WASP-South 
lightcurves show no indication of any ellipsoidal 
variation due to the distortion of the star by a 
massive companion. 

We obtained 22 radial velocity measurements 
during the interval 2010 January 3 to 2010 August 
5 with the fibre-feb CORALIE spectrograph on the 
Euler 1.2-m telescope located at La Silla, Chile. 
The spectra were obtained with an exposure time 
of 30 minutes and have a typical signal-to-noise 
ratio of 25-30. Accurate wavelength calibration 
is ensured by the simultaneous observation of a 
thorium-argon arc in a second fibre feed to the 
spectrograph. Details of the instrument and data 
reduction can be found in lQueloz et al. and 
references therein. The RV measurements were 
performed using cross-correlation against a nu- 
merical mask generated from a G2-type star are 
given in Table [H where we also provide the bisec- 
tor span, BS, which measure s the asymmetry o f 
the cross-correlation function (jQueloz et al. II2OOII) . 
The standard error of the the bisector span mea- 
surements is 2crRv- 

We also obtained photometry of TYC 7247-587- 
1 and other nearby stars on 2010 June 23 using 
the LCOGT 2.0-m Faulkes Telescope South (FTS) 
at Siding Spring Observatory. The Merope cam- 
era we used has an image scale of 0.279 arcsec- 
onds/pixel when used in the 2x2 binning mode we 
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employed. We used a Pan- STARR^I z-band filter 
to obtain 210 images covering one transit. These 
images were processed in the standard way with 
IRAlil using a stacked bias image, dark frame, 
and sky flat. T he DAOPHOT photometry package 
(jStetson 1987 ) was used to perform object detec- 
tion and aperture photometry for WASP-41 and 
several comparison stars in the 5' x 5' field of view 
of the instrument. Observations were interrupted 
by poor weather so there are no observations dur- 
ing the egress phase of the transit. These data are 
sufficient to confirm that the transit-like features 
seen in the WASP-South data are due to the star 
TYC 7247-587-1 and to provide better measure- 
ments of the depth of the transit than is possible 
from the WASP-South data (Fig.©. 

All photometric data presented in this paper 
are available from the NStED databaseH 

3. WASP-41 Stellar Parameters 
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Fig. 3. — Section of the co-added CORALIE spectra 
of WASP-41 showing emission in the core of the 
Call H and K lines. 



The 17 individual CORALIE spectra of WASP- 
41 available up to May 2010 were co-added to pro- 
duce a single spectrum with a signal-to-noise ratio 
»70:1. The standard pipeline reduction products 
were used in the analysis. 

The a nalysis was p e rform ed using the methods 



given m 



Gillon et al.l (|2009[) . The line was 
used to determine the effective temperature (Toff), 
while the Na i D and Mg i b lines were used as sur- 
face gravity (log 5) diagnostics. The parameters 
obtained from the analysis are listed in Table [5] 
The elemental abundances were determined from 
equivalent width measurements of several clean 
and unblended lines. A value for micro-turbulence 
H i) was d e termi ned from Fe i using the method 
of lMagaiii ( 19841 ). The quoted error estimates in- 



clude that given by the uncertainties in Tcs, log 5 
and ^t, as well as the scatter due to measurement 
and atomic data uncertainties. 

The projected stellar rotation velocity (usinz) 
was determined by fitting the profiles of several 
unblended Fe i lines. A value for macro-turbulence 
(I'mac) of 2.3 ± 0.3 kms~^ was assumed, based on 
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Fig. 5. — Radial velocity and bisector span mea- 
surements for WASP-41. Upper panel: Radial ve- 

^http: //pan- Starrs ■ if a.hawail . edu/public/deslgn-f eatures/cam 4g^tKtiH?lative to the centre-of-mass velocity to- 
^IRAF is distributed by the National Optical Astronomy gether with the bcst-fit circular Orbit. Lower 

Observatory, which is operated by the Association of Uni- panel: bisector span measurements, 

versities for Research in Astronomy (AURA) under coop- 
erative agreement with the National Science Foundation, 
• ^http: //nsted. Ipac . caltech. edu| 
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Fig. 4.— Periodogram of the WASP data from two seasons for WASP-41 and TYC 7247-459-1. The star 
name and date range (JD— 245000) are given in the title of each panel. Horizontal lines indicate false alarm 
probability levels FAP=0. 1,0. 01, 0.001. Note the change of scale in the lower right panel. 



the tabulation bv lGrav ( 20081 ). and an instrumen- 
tal FWHM of 0.11 ± 0.01 A, determined from the 
telluric lines around 6300A. A best fitting value of 
ijsini — 1.6 ± 1.1 kms^^ was obtained. 

We observe emission lines in the cores of 
the Call H and K lines in our spectrum of 
WASP-41 (Fig. |3|). We estimate a s pectral index 
log R '„ y « —4.67 from this spectrum (jNoves et al 



19841 ) , but note that a transformation for measure- 
ments made with the CORALIE spectrograph to a 
standard system does not yet exist. 

4. Rotation period 

The presence of Call H and K emission lines 
in the spectrum of WASP-41 suggests the possi- 
bility of variations in the lightcurve due to star 
spots with the same period as the rotation period 
of the star. The projected equatorial rotation ve- 
locity of WASP-41 combined with the estimated 
radius given in Table [2] imply a rotation period of 
32 ± 23 d. We therefore searched for photometric 
variations in our WASP photometry with a similar 
period. 



Our WASP photometry has a large number of 
observations (several thousand) with standard er- 
rors that vary from a few milli-magnitudes for 
the best data up to a magnitude for data ob- 
tained on cloudy nights. WASP data obtained 
on cloudy nights is sometimes affected by large 
systematic errors, so we remove data with stan- 
dard errors larger than 5 times the median value 
prior to further analysis. The remaining data still 
have a range of standard errors to it is advanta- 
geous to use a period searching algorithm that can 
account for this. The generalized Lomb-Scargle 
periodogram defined by IZechmeister and Kiirster 
(12009,1 is a suitable method. However, there is lit- 
tle to be gained from including a "floating mean" 
in the case of the WASP data, so we used a slightly 
different definition of the periodogram equivalent 
to a least-squares fit of the sinusoidal function 
i/i = asin(cLiti) -|- bcos{ujti) to magnitudes = 
1, 2, . . . , with with standard errors Ci obtained 
at times U for a given angular frequency w. Note 
that we implicitly assume that the magnitudes 
have a weighted mean value of 0. We used the fol- 
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Table 1: Radial velocity measurements for WASP- 
41 



BJD 




RV 




crv 




RS 




-2 450 000 


(km 


s-i) 


(km 




(km 


s-1) 


5200 


8101 


3 


1691 





0076 


-0 


0168 


5290 


8090 


3 


3837 





0070 


-0 


0323 


5293 


8035 


3 


4106 





0087 


-0 


0369 


5294 


8225 


3 


2084 





0077 


-0 


0193 


5295 


7751 


3 


1732 





0069 


-0 


0048 


5296 


7550 


3 


3926 





0070 


-0 


0240 


5297 


8185 


3 


2263 





0083 


-0 


0076 


5298 


7732 


3 


1696 





0073 


-0 


0303 


5300 


8114 


3 


2652 





0084 


-0 


0067 


5301 


8036 


3 


1611 





0088 


+0 


0268 


5305 


7269 


3 


3935 





0071 


-0 


0239 


5306 


7188 


3 


3411 





0072 


-0 


0043 


5311 


8040 


3 


3486 





0088 


-0 


0356 


5317 


5927 


3 


3310 





0101 


-0 


0144 


5323 


6579 


3 


2985 





0073 


-0 


0193 


5326 


6585 


3 


2569 





0071 


+0 


0065 


5327 


6832 


3 


4099 





0133 


-0 


0278 


5362 


4677 


3 


1436 





0123 


+0 


0065 


5375 


5148 


3 


2806 





0080 


-0 


0043 


5388 


5621 


3 


4714 





0092 


-0 


0490 


5410 


5201 


3 


3512 





0083 


-0 


0182 


5414 


4912 


3 


1813 





0083 


-0 


0024 



Table 2: Stellar parameters of WASP-41 from 
spectroscopic analysis. 



Parameter 


Value 


Jeff [J^-J 


040U ± iOU 




4.4 it U.Z 


?t [kms J 


l.U ± U.z 


vsiiii [kms~^] 


i.D ± i.i 


[i'e/HJ 


n HQ 4- n no 
— u.Uo ± u.uy 


\~W^ /TT1 

[i\ a/ rlj 


n n7 -1- A nn 
U.U/ ± U.U9 


[Mg/MJ 


U.U / ± U. ID 


r A 1 /xji 
[Al/HJ 


A A1 _L A AO 

— U.Ui ± U.Uo 


[bi/HJ 


A AC _L A t\Ci 

U.Uo ± U.Uo 


fPo /Wl 


U.Uo ± U.iO 


fcip /"Rl 
[oc / ni 


A ni -1- n 1 

u.ui in u.iu 


[Ti/H] 


0.00 ± 0.10 


[V/H] 


0.06 ± 0.17 


[Cr/H] 


0.00 ± 0.05 


[Mn/H] 


0.00 ± 0.13 


[Co/H] 


-0.01 ± 0.07 


[Ni/H] 


-0.04 ± 0.06 


log A(Li) 


< 0.5 


Mass [Mq] 


0.95 ± 0.09 


Radius [Rq] 


1.01 ± 0.26 


Sp. Type 


G8 V 


Distance [pc] 


180 ± 60 



Note: Ma. s s and Radius estimate using the 
Torres et al. ( 2010t ) calibration. Spe ctral Type e s- 
timated from Toff using the table in iGravl (|2008l ). 



lowing definition of the power at angular frequency 



xl ' 



where 



and 



2 sr^^i 

(wj - y,)2 



We used this definition of the power to search 
for periodicity in the WASP-South lightcurves of 
WASP-41 at 4096 evenly spaced frequency val- 
ues from to 2.5 cycles/ day. Variability due to 
star spots is not expected to be coherent on long 
timescales as a consequence of the finite lifetime 
of star-spots and differential rotation in the pho- 
tosphere so we analysed the two seasons of data 
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for WASP-41 separately. The signal of the tran- 
sits in the data was removed using a model similar 
to the one described below, although in practice 
this has a negligible effect on the periodogram. We 
also analysed 3 nearby stars of similar magnitude 
and colour observed simultaneously with the same 
camera. The results are shown in Table [S The 
periodograms of WASP-41 and one of the nearby 
stars are shown in Fig.H 

We used the method of IPress and Rvbicki 



()l989l) ;o allow for fast and accurate computa- 
tion of the periodogram. This enabled us to 
investigate the significance levels of any peaks 
in our periodograms using a boot-strap Monte- 
Carlo approach. We generated synthetic data 
sets using a method s i milar to that devised by 



Collier Cameron et al 



(l2009l) . in which data are 
shuffled according to the date on which the data 
were obtained. The shuffiing procedure shifts each 
night's observations in their entirety to a new date. 
This procedure effectively destroys coherent sig- 
nals with periods longer than 1 day, but retains 
the global form of the window function and the ef- 
fects of correlated noise. In additional to shuffling 
the date of the observations, we also change the 
sign of all the data from a given night for a random 
selection of half of the dates of observation. 



Collier Cameron et al. ( 20091 ) use a single sim 



ulated data set to estimate the parameter N^s in 
the following model for the distribution of power 
in the periodogram in the absence of any periodic 
signal; 

Prob(P„>P;) = (l-i^;)(^--^)/^ 

They then use this model to estimate the false 
alarm probability, FAP, for the highest peak in 
the periodogram of the actual data from 



FAP = 1 - [1 - Prob(P„ > F^)] 



where M ~ {tpj — ti)(/max — /min) is an estimate 
of the number of independent frequencies in the 
periodogram calculated over the range of frequen- 
cies /mill to /maxH We found from simulations 
of a large number of data sets using both the 
"night-shuffling" method and Gaussian random 



^Note that there is an error in equation (2) of 
ICoIIier Cameron et al.L the equation for FAP given here is 
the one actually used in that paper. 



noise that this method gives inaccurate estimates 
of FAP, particularly when this value is small. The 
principle reason for this is that the value of Noff 
estimated from a single periodogram has a rather 
large uncertainty, typically around 20% for the 
data sets we investigated. We found that the fol- 
lowing empirical model gives a good representa- 
tion of the distribution of the peak power Pn.bcst 
in the periodograms we generated from simulated 
data 



= Co 



ciP' 



best 



+C2 ( ^'nhost 



logio[Prob(P„,bcst > ,bcst / 



For the analysis of a given set of data we use this 
model to fit the distribution of peak power val- 
ues in 1024 simulated data sets. We then use the 
values of cq, ci and C2 to calculate the power cor- 
responding to FAP=0.1, 0.01 and 0.001 shown in 
Fig. U and to estimate the FAP value for the high- 
est peak in the actual periodogram given in Ta- 
ble El 

There is a clear peak in the periodogram for 
WASP-41 near P=18d, although there appears to 
be a difference in the period derived from the two 
seasons of data (17.6 d and 18.4 d). We inspected 
the periodograms from the first season of data for 
WASP-41 and the three nearby stars and found 
that they all show power near Id, 30 d and var- 
ious combinations of these frequencies and their 
harmonics, presumably as a result of systematic 
errors in the photometry related to the diurnal 
and lunar cycles. In contrast, there is very little 
spurious power in the periodograms for the sec- 
ond season of data. For this reason we identify 
the period of 18.41±0.05d derived from the sec- 
ond season of data as the correct rotation period 
for WASP-41. We estimated the standard error of 
the period measurement by analysing 1024 data 
sets with the same number of points as the origi- 
nal sample, randomly re-sampled with reselection 
from the original data. 

5. Planetary parameters 

Our radial velocity (RV) and bisector span (BS) 
measurements for WASP-41 are shown in Fig. [5] 
as a function of the photometric transit phase. 
There is a weak but statistically significant anti- 
correlation between these RV and BS values. In 
principle, such an anti-correlation can be used to 
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Table 3: Period analysis of WASP-42 and nearby stars of similar magnitude and colour. 



Star 


V 


B-V 


Year 


N 


RMS 


Pbcst/rf 


Power 


FAP 


WASP-41 


11.6 


0.7 


2007 


5673 


0.041 


17.62 


0.358 


2 X 10-5 








2008 


4966 


0.019 


18.41 


0.240 


6 X lO"'^ 


TYC 7247-1008-1 


11.8 


0.3 


2007 


5414 


0.026 


1.00 


0.167 


0.07 








2008 


4834 


0.012 


0.98 


0.023 


0.65 


TYC 7247-459-1 


11.5 


0.5 


2007 


5389 


0.025 


96.38 


0.101 


0.07 








2008 


5003 


0.012 


1.62 


0.006 


0.45 


TYC 7247-683-1 


11.2 


0.5 


2007 


5325 


0.024 


86.23 


0.089 


0.06 








2008 


4993 


0.012 


0.46 


0.005 


1.00 



identify spurious RV signals caused by st ellar ac- 
tivity. For example, (jQueloz et al.l 120011) showed 
that the apparent RV signal in HD166435 is caused 
by stellar activity by identifying 3 characteristics 
of the signal: i. the RV signal was not coher- 
ent over timescales of more than 30 days; ii. the 
amplitudes of the BS and RV variation were sim- 
ilar; iii. the BS variations were correlated with 
the photometric variations of the star. In the case 
of WASP-41 there is no doubt that the RV signal 
is due to the presence of a planetary mass com- 
panion to WASP-41 and not due to stellar activ- 
ity because: i. the RV signal is coherent over a 
timescale of more than 200 days; ii. the ampli- 
tude of the BS variation is an order of magnitude 
smaller than the RV variation; iii. the BS values 
are also correlated with the phase calculated for a 
rotation period of 18. 4d established from the pho- 
tometric variations of WASP-41. As an additional 
test of the reality of the planetary RV signal we 
also re-measured the RV values using a K5-type 
mask instead of a G2-type mask. Looking for a 
difference between the RV values derived with the 
different masks is an effective method for detecting 
spurious RV signals caused by star spots or faint 
eclipsing binary stars blended with the target star. 
In the case of WASP-41 there is no significant dif- 
ference in the RV values measured with the two 
masks so the best explanation for the observed RV 
and BS signals is that WASP-41 has a planetary 
mass companion with an orbital period of 3.05 d. 

The CORALIE radial velocity measurements 
were combined with the WASP-South and FTS 
photometry in a simultaneous Markov-chain Monte- 
Carlo (MCMC) analysis to find the parameters 
of the WASP-41 system. We removed the pho- 
tometric variation due to stellar rotation from 



the WASP-South photometry by subtracting a 
first-order harmonic series, fit by least-squares 
to the differential magnitudes from each season 
independently. The shape of the transit is not 
well defined in the WASP-South or FTS pho- 
tometry, so we have imposed an assumed main- 
sequence mass-radius relation as an additional 
constraint in our analysis of the data. The stel- 
lar mass is determined from the parameters Tcs, 
logq and [Fe / H] us ing the procedure described by 
Enoch et al.l (l2010t ). bas ed on the cq i npilat ion of 



eclipsing binary data bv lTorres et al. I (|2010l) . The 
code uses T^s and [Fe/H] as MCMC jump vari- 
ables, constrained by Bayesian priors based on the 
spectroscopically-determined values given in Ta- 
blej2l Limb-darkening coefficients are taken from 
IClaretl (l2Q00l) . 



The parameters derived from our MCMC anal- 
ysis assuming a circular orbit are listed in Tabled 
We found that the contribution to the value of 
for the fit from the RV data was much higher than 
expected given the standard errors of these mea- 
surements. This is a result of additional noise in 
the RV data due to stellar activity ("jitter"). The 
free parameters in the fit to the RV are Tq, P, 
7 and K, but the values of Tq and P are deter- 
mined almost entirely by the photometric data, so 
the number of degrees of freedom in the fit to the 
RV data is approximately 20. To achieve a con- 
tribution to the value of from our 22 RVs of 
approximately 20 we found that we needed to add 
21 ms"^ in quadrature to the standard errors for 
the RV data. We also performed an MCMC anal- 
ysis of the data including the parameters -^/ecos w 
and -yesinw as free parameters. These parame- 
ters are used to describe an eccentric orbit solu- 
tion because they are not strongly correlated and 
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are equivalent to a assuming uniform prior distri- 
bution for the value of the eccentricity, e. With 
this eccentric-orbit solution and the same value 
of the jitter as the circular orbit solution we find 
e = 0.10 ± 0.06. As this is consistent with a cir- 
cular orbit we adopt the parameters from the cir- 
cular orbit solution. As we do not know a priori 
that the orbit is circular, we take the standard er- 
rors on the parameters from the non-circular orbit 
solution. 

We considered the contribution of correlated er- 
rors ("red noise") to the standard errors quoted 
in Table ID While the individual transits in the 
WASP data are affected by red noise, the anal- 
ysis of the combined lightcurve covering many 
individual transits will not be strongly affected 
by red noise because there will be no correla- 
tion between the systematic noise from different 
nights. The FTS lightcurve is affected by red noise 
so we have investigated the effect of this using 
the "prayer bead" method. A separate MCMC 
analysis was performed in which synthetic FTS 
lightcurves were created from the model fit to the 
lightcurve and the residuals from this model af- 
ter cyclic permutation at each step in the MCMC 
chain. We find that this does not significantly in- 
crease the error estimates for any of the parame- 
ters. 

Another issue to consider for the FTS lightcurve 
is the effect of starspots on the lightcurve, par- 
ticularly given the sparse phase coverage of this 
lightcurve. Star-spots can affect the parameters 
derived from this lightcurve in two ways. Firstly, 
the normalization of the lightcurve can be affected 
by the overall change in brightness of the system 
due to the 18.4 d rotational modulation of the 
lightcurve. Secondly, small star spots covered by 
the planet during the transit by distort the shape 
of the lightcurve in a way that may be hard to 
spot by-eye, but that introduc e a systematic er- 



ror in the parameters derived (jMiller-Ricci et al 
20081 ). We have identified parameters in our so- 



lution that are determined primarily by the FTS 
lightcurve by performing an MCMC solution us- 
ing the WASP and CORALIE data only. We find 
that the depth of transit and the planet radius are 
the only parameters for which the standard er- 
ror estimates increase significantly in the solution 
excluding the FTS data. In both cases, the val- 
ues from the two solutions agree, but the solution 



excluding the FTS data has standard error esti- 
mates approximately twice as large as the solution 
including the FTS data. We conclude that these 
two parameters may be subject to systematic er- 
rors comparable to the random errors quoted in 
Table H 

We have compared the values of the stellar ef- 
fective temperature, T^g, a nd the stellar density , 

(l2000h . 



to the stellar models of Girardi et al. 



P 

We use the parameters Tf,g and because they 
are independently determined directly from the 
observations. We find that the mass inferred from 
the models (0.9 ± O.IM0) is consistent with the 
mass derived in our MCMC analysis and that the 
uncertainties on the values of Toff and are too 
large for the models to provide any useful con- 
straint on age of the star. The surface gravity de- 
rived from our MCMC solution is consistent with 
the log g value from the analysis of the spectrum, 
but the large uncertainty on the latter value means 
that this is a rather weak constraint. 

6. Discussion 

WASP-41b joins a growing number of plan- 
ets discovered with masses ss 0.9Mjup, radii w 
1.2i?jup orbiting solar- like stars with periods of 
about 3 days, e.g., HAT-P-13b, XO-lb, WASP- 
28 b, CoRoT-12b, WASP-26b, HAT-P-5b, HAT- 
P-6b, etcH The WASP-41 planetary system is 
also similar to the TrES-1 system, particularly in 
regard to the activity level of the host star. 

Knutson et al. ( 2010l) suggest that there is a 



clear connection between the properties of the 
emission spectra from hot Jupiters and the ac- 
tivity levels of their host stars. In general, the 
emission spectrum of hot Jupiters measured from 
the eclipse depths at infrared wavelengths can 
only be matched by models that include a high- 
altitude temperature inversion. Four exceptions 
to this general rule are the planets orbiting the 
stars HP 1897 33, TrES-3, TrES-1 and WASP-4. 
Knutson et al. found that all four of these stars 



show moderate levels of chromospheric activity, 
higher than for all the other stars for which they 
were able to measure a value log R'lj^j^ ■ The avail- 
able evidence supports a hypothesis in which some 
source of optical opacity high in the atmosphere 



^ http:/ / exoplanet. eu\ 
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of hot Jupiter planets causes a temperature inver- 
sion, but that this opacity source is destroyed by 
the UV flux associated with chromosph eric activ- 
ity. The result found bv lKnutson et al. has a high 
statistical significance but is based on a sample of 
only 15 stars. The situation is further complicated 
by the correlation between a planet's surface grav- 
ity and the level of c hromosph e ric ac tivity in its 
host star claimed bv iHartman ( 201Clf l. The val- 
ues of log R'uii and log Qp we have measured for 
WASP-41 are consistent with the correlation ob- 



served bv lHartmanl . The connection between chro- 
mospheric activity and temperature inversions can 
be confirmed in the case of WASP-41 using obser- 
vations of the secondary eclipse with the IRAC in- 
strument o n the Spitzer Space Te lescope at 3.6/im 
and 4.5/im (jKnutson et al.ll2QlCll ). 

In principle, one can use an age -activity re- 
lation established from the average behaviour of 
many solar-type stars to estimate the age of a star 
based on the measured log R'^}^ value. However, 
solar-type stars have activity cycles with periods of 
about a decade during whic h the value o f log R'ut( 
can vary by 10% or more (|Baliunas et al. 1995). 
This means that an age estimate based on a sin- 
gle measurement of log R'^j^ has an unknown sys- 
tematic error that can be large enough to make 
the age estimate effectively meaningless. A more 
useful age estimate can be made in this case 
based on the observed rotation period of the star. 
Plot — 18.4 d. The c alibration of th e "gyrochrono- 
logical" age given by ^ Barnes ( 2007 ) implies an age 
of 1.8 Gyr for WASP-41 with an error of about 
15%. Stars with Toff =5300- 5600K are seen to be 
lithium poor in clusters with ages > 1 Gyr such as 
NGC752, M67 but the majority of stars in this Tcff 
range in clusters with ages of about 600 Myr such 
as Praesepe and the Hyades, NGC6633 and Coma 
Berenices are lithium rich. We conclude that the 
upper limit to the lithium abundance given in Ta- 
ble [5] is consi stent with any age greater than about 
600 Myr (,Sestito and Randich,2005.) . 

WASP-South is hosted by the South African 
Astronomical Observatory and we are grateful for 
their ongoing support and assistance. Funding 
for WASP comes from consortium universities and 
from the UK's Science and Technology Facilities 
Council. We thank the referee for their careful 
reading of the manuscript and for comments that 



improved the quality of this paper. 
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Table 4: System parameters for WASP-41. The planet equilibrium temperature is calculated assuming a 
value for the Bond albedo A= 0. N.B. an assumed main-sequence mass-radius relation is imposed as an 
ad ditional const r aint in this solution so the mass and radius of the star are not independent parameters - 
Enoch et al. I (l201(j) for details. Parameter values are taken for the solution assuming e = 0, standard 



see 



errors on the parameters are taken from the solution with y^cos(a;) and ■\/esin(a;) as free parameters. 



Parameter 


Symbol 


Value 


Units 


Transit epoch (HJD) 


To 


2455343.463 ±0.001 


days 


Orbital period 


P 


3.052401 ± 0.000004 


days 


Planet/star area ratio" 


[Rpl Ri,)"^ 


0.0186 ±0.0004 




Transit duration 


tr 


0.108 ±0.002 


days 


Impact parameter 


b 




R^ 


Stellar reflex velocity 




0.135 ±0.008 


kms^^ 


Centre-of-mass velocity 


7 


3.284 ±0.009 


kms^^ 


Orbital eccentricity 


e 


(fixed) 




Orbital inclination 


i 


87.7 ±0.08 


degrees 


Stellar density 


P* 


1.27±0.14 


Pq 


Stellar mass 




0.93 ±0.03 




Stellar radius 


i?* 


0.90 ±0.05 


Rq 


Orbital semi-major axis 


a 


0.0402 ± 0.0005 


AU 


Planet radius" 


Rp 


1.20 ±0.06 


Rj 


Planet mass 


Mp 


0.92 ±0.06 


Mj 


Planet surface gravity 




3.16 ±0.04 


[cgs] 


Planet density 


Pp 


0.50 ±0.08 


Pj 


Planet temperature 


T 

^ cq 


1230 ±50 


K 



" These quantities may have a systematic error comparable to the random error due to the influence of stars 
spots on the FTS lightcurve. 
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This 2-column preprint was prepared with the AAS lAlJ^X 
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